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Abstract 

If ultra-high energy cosmic rays (UHECRs) are accelerated at astrophysical point 
sources, the identification of such sources can be achieved if there is some kind of 
radiation at observable wavelengths that may be associated with the acceleration 
and/or propagation processes. No radiation of this type has so far been detected or 
at least no such connection has been claimed. The process of photopion production 
during the propagation of UHECRs from the sources to the Earth results in the 
generation of charged and neutral pions. The neutral (charged) pions in turn decay 
to gamma quanta and electrons that initiate an electromagnetic cascade in the 
universal photon background. We calculate the flux of this gamma radiation in 
the GeV-TeV energy range and find that for source luminosities compatible with 
those expected from small scale anisotropics in the directions of arrival of UHECRs, 
the fluxes can be detectable by future Cerenkov gamma ray telescopes, such as 
VERITAS and HESS, provided the intergalactic magnetic field is not larger than ~ 
10^^*^ Gauss and for source distances comparable with the loss length for photopion 
production. 
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1 Introduction 



The discovery of the cosmic microwave background radiation (CMBR) [1] has 
been a major breakthrough not only for cosmology but for cosmic ray physics 
as well. Soon after its discovery, it was shown that cosmic rays would suffer 
severe energy losses due to the inelastic process of photopion production, pro- 
vided the energy of cosmic rays being higher than the kinematic threshold 
for this process. This effect was predicted to result in a flux suppression at 
energies around 10^° eV [2,3] for a homogeneous distribution of the sources, 
something which is now known as the GZK feature. Despite the efforts ded- 
icated to the experimental detection of this feature, to date we do not have 
a definite answer to whether the GZK suppression is there or not. The two 
largest detectors currently operating, AGASA and HiRes, give results which 
appear discrepant, although the statistics of events does not allow us to infer a 
definitive conclusion on this issue [4]. The detection of the GZK feature in the 
spectrum of UHECRs would be the clearest proof that UHECRs are generated 
at extragalactic sources. 

Despite the efforts in building several large experiments for the detection of 
particles with the highest energies in the cosmic ray spectrum, the nature 
of the sources of these particles and the acceleration processes at work are 
still unknown. No trivial counterpart has been identified by any of the current 
experiments: there is no significant association with any large scale local struc- 
ture nor with single sources. The lack of counterparts is particularly puzzling 
at the highest energies, where the loss length of the detected particles is small 
enough that only a few sources can be expected to be located within the error 
box of current experiments. 

It has been recently proposed that a statistically significant correlation may 
exist between the arrival directions of UHECRs with energy above 2 x 10^^ eV 
and the spatial location of BL Lac objects, with redshifts larger than 0.1 [5,6] 
(but see also [7] for a criticism and [8] for the reply). Since at these energies 
the loss length is comparable with the size of the universe, no New Physics 
would be needed to explain the possible correlation. On the other hand, since 
no BL Lac is known to be located close to the Earth, the spectrum of UHECRs 
expected from these sources would have a very pronounced GZK cutoff. 

The potential for discovery of the sources has recently improved after the iden- 
tification of a few doublets and triplets of events clustered on angular scales 
comparable with the experimental angular resolution of AGASA. A recent 
analysis of the combined results of most UHECR experiments [9] revealed 8 
doublets and two triplets on a total of 92 events above 4 x 10^^ eV (47 of which 
are from AGASA). The recent HiRes data do not show evidence for significant 
small angle clustering, but this might well be the consequence of the smaller 
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statistics of events and of the energy dependent acceptance of this experiment: 
one should remember that in order to reconstruct the spectrum of cosmic rays, 
namely to account for the unobserved events, a substantial correction for this 
energy dependence need to be carried out in HiRes (the acceptance is instead a 
flat function of energy for AGASA at energies above 10^^ eV) . This correction 
however does not give any information on the distribution of arrival directions 
of the missed events. 

The statistical significance of these small scale multiplets of events have re- 
cently been questioned in [10]. If the appearance of these multiplets in the data 
will be confirmed by future experiments as not just the result of a statistical 
fluctuation or focusing in the galactic magnetic field [11], then the only way 
to explain their appearance is by assuming that the sources of UHECRs are 
in fact point sources. This would represent the first true indication in favor 
of astrophysical sources of UHECRs, since the clustering of events in most 
top-down scenarios for UHECRs seems unlikely. A clear identification of the 
GZK feature in the spectrum of UHECRs would make the evidence in favor 
of astrophysical sources even stronger. 

In some recent work [12,13,14,15,16] it has been shown that the small scale 
anisotropies can be used as a powerful tool to measure the number density 
of sources of UHECRs. In particular, in [12] the number density of sources 
that best fits the observed number of doublets and triplets was found to be of 
order 10"^ - 10"^ Mpc'^ with an average luminosity per source of 2 x 10 — 
2 X lO^^erg s~^ at energies in excess of 10^^ eV. The extrapolation of these 
numbers down to GeV energies are strongly dependent upon the injection 
spectrum, assumed here in such a way to fit the observed spectrum of cosmic 
rays above 10^^ eV. 

Alternative estimates of the density of sources of UHECRs recently appeared 
in the literature: the number obtained in [17] is in fair agreement with the re- 
sult of [12], while more complex is the comparison with the case investigated 
in [18,46]. In [18], a source density of ~ 10^^ — 10~^Mpc~'^ was derived, while 
the large scale isotropy of UHECRs was used to infer a local field strength of 
~ 0.1//G in the local supercluster. This conclusion was mainly due to the fact 
that the authors neglected the contribution of distant sources of cosmic rays 
with energy above 4 x 10^^ eV. When this effect is accounted for, as in a more 
recent paper [46] by the same authors, the revised estimate of the source den- 
sity appears to be in rough agreement with that of [12] for the cases of weak 
magnetization. The authors of [46] also investigate the cosmic ray propagation 
in a magnetized universe, as obtained from numerical simulations of large scale 
structure formation with passive magnetic fields. In these cases the estimate 
of the source density becomes more model dependent. It is probably worth 
stressing however, that the magnetic fields obtained in [46] do not appear to 
be compatible with those derived in the numerical MHD simulations of [45]: 
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while the magnetic field obtained inside clusters of galaxies in both numerical 
approaches appears to have roughly the same strength, the fields outside clus- 
ters differ by 2-4 orders of magnitude. Not surprisingly the conclusions that 
the two groups draw on the propagation of UHECRs are incompatible with 
each other. It would be auspicable that the two pieces of work can achieve 
some kind of agreement, so to clarify the situation. 

The interactions of cosmic ray protons at the highest energies with the photons 
in the CMBR, responsible for the appearance of the GZK feature in the dif- 
fuse spectrum of UHECRs, also generate gamma radiation through the decay 
of neutral pions and electrons through the decay of charged pions (actually 
mostly positrons are produced, but hereafter we will refer to both electrons 
and positrons as electrons). The very high energy gamma radiation produced 
in this way initiates an electromagnetic cascade that results in the appearance 
of a gamma ray fiux in the energy region accessible to Cerenkov detectors. 
We investigate in detail this process and assess the detectability of this signal 
from the directions where the sources of UHECRs are located. 

The paper is organized as follows: in section 2 we present an analytical estimate 
of the development of the electromagnetic cascade initiated by UHECRs; in 
section 3 we detail the background radiation field we used for our calculations; 
in section 4 we present the detailed calculation of the electromagnetic cascade 
developement; in section 5 we present our results for different luminosities of 
the sources and different values and topologies of the extragalactic magnetic 
field; in section 6 we draw our conclusions. 



2 The electromagnetic cascade initiated by UHECRs: an analytical 
approach 

The purpose of this section is to provide an estimate of the expected fiux of 
gamma rays due to the cascade initiated by UHECRs. A detailed calculation 
will be presented in the next sections, where many effects neglected here will 
be properly accounted for. 

UHECRs with energy above the threshold for photopion production are kinc- 
matically allowed to interact inclastically with the photons in the CMBR and 
generate charged and neutral pions mainly through the reactions 

p -|- 7 — > TT^ -|- n p + ^^n^+p. (1) 

The decay of neutral pions results almost instantaneously in the production 
of two photons with typical energy ~ ~ O.lEp (here, for the purpose 
of a rough estimate we assumed an inelasticity of ~ 20%). 
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These photons start an electromagnetic cascade through pair production on 
the universal photon background. A similar cascade is initiated by positrons 
resulting from the decays of tt+'s. 

Cosmic rays with energy above ~ 4 x 10^^ eV can interact with the CMBR 
through proton pair-production, which results in electrons and positrons which 
can also initiate an electromagnetic cascade. This process becomes however 
relevant for the propagation of UHECRs only at energies above ~ 10^^ eV, 
when adiabatic energy losses due to the expansion of the universe become 
less severe. We estimated that the contribution of the proton pair production 
to the electromagnetic cascade is of few percent of the contribution initiated 
through photopion production but can be up to ~ 50% for distant sources. In 
the rest of this calculation wc shall focus on the contribution from photopion 
production, while a full treatment of additional contributions will be properly 
considered in an upcoming paper (Ferrigno, Blasi & De Marco, in preparation). 

We will provide a detailed description of our calculations of the development 
of the cascade in the next section. Here we estimate the expected effect using 
an approach which was proposed in [19] and detailed in [20]. 

Let us consider a simple situation in which the cosmic photon background is 
dominated by the radio radiation and the CMBR. For our purposes here we 
can assume that the typical energy of the radio photons is eradio = 4 x 10~^ 
eV (corresponding to a frequency of 1 MHz). For the CMBR, the reference 
energy will be taken as ecus = 6 x 10~^ eV. 

Most of the gamma rays that initiate the electromagnetic cascade have energy 
E-y > 10^^ eV. We define here the two quantities Xradio — Eeradio/nT'l and 
xcMB — EecMB l^li where E is the energy of the leading particle, which in 
our case starts as being either a photon or an electron. Clearly Xradio ^ xqmb 
therefore the cascade develops first on the radio photon background. For the 
particles that start the cascade, Xradio ~ 0.1 — 1. In fact Xradio can be quite 
higher for higher energy photons, that are also generated as a result of the 
photopion production. For the range of values Xradio ^ 1 the pair production 
of photons and the inverse Compton scattering (ICS) of electrons once they 
have been produced, occur in the Klein-Nishina (KN) limit and the leading 
particle loses a fraction of its energy which is approximately / = AE/E pa 
1 / l-a(2xradio) ■ The energy of the secondary particle produced will therefore be 
E' fE. The KN regime stops however when Xradio ~ 1, which is almost 
immediately the case for the interaction on the radio background. At this 
point the cascade develops efficiently through pair production and ICS down 
to the point when photon energies fall below the threshold for pair production, 
^th,radto _ ml/ eradio, and the number of electrons remains constant at later 
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times. Therefore we can write for the electron spectrum 



qe{Ee) = go if Ee < E] 



ith, radio 
'7 



(2) 



At higher energies, conservation of energy gives an electron spectrum that 
declines as the first power of the electron energy [19,20]. The photon spectrum 
in the cascade is generated through the ICS of these electrons. Since the ICS 
occurs in the Thomson regime, it is simple to demonstrate that the electrons 
with a fiat spectrum (ge(-E'e) = qo) generate a gamma ray spectrum n-y{E^) oc 
E'^/"^ up to a maximum energy 



The electrons with spectrum qe{Ee) oc E~^ generate inverse compton photons 
with spectrum n^{E^) oc E'"^ up to the energy = rnl/eradio- At 

higher energies there are no gamma rays because the electron spectrum is 
cut off. The normalization in the gamma ray spectrum is such that the total 
energy in it equals the energy in the primary particles (gamma rays and/or 
electrons) . 

It is worth noticing that the spectrum of the initial gamma rays and electrons 
does not explicitely enter the calculation: the spectrum of the final cascade 
is universal. In other words the electromagnetic cascade behaves as a sort of 
calorimeter that redistributes the initial energy into gamma rays with a given 
spectrum. Violations to this rule will be seen and discribed in the next section, 
but the basic result is not substantially changed. 

Once this first part of the cascade has developed, the gamma rays generated 
so far can produce a second cascade on the CMB photon background. This 
second cascade follows exactly the same steps described above. We introduce 
here a new energy scale 



The electron spectrum is flat up to the energy £'*^.ca^^/2 as in the first part of 
the cascade. The parent gamma ray spectrum, as we found above, has the first 
change in slope at the energy E'x'^'^", therefore it follows that a slope change 
will appear in the electron spectrum generated through pair production, at 
the energy ^^^"^^72. At energies E^J^'^^^^ /2 < E^ < E''^'^''"/2 the electron 
spectrum is qe{Ee) oc E~^^'^, while at higher energies it again steepens to 
qe{Ee) OC E~^. The final gamma ray spectrum can be calculated as ICS of the 




(3) 
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electrons, and after simple calculations we obtain: 




rp ^ rpCMB 



A — 7 — 3 mj. 



(f) 



2 



E<E^< E'J"'^^^ 



(5) 



The constant no should be calculated by requiring that the total energy of 
the primary gamma rays and electrons is converted into the energy of the 
cascade. The calculations outlined above assume in principle that we start 
with a gamma ray and/or electron beam at some distance. In fact, for the 
case of gamma rays generated from the decay of neutral pions created in 
turn in the p + 7 inelastic interactions this is not the case: gamma rays and 
electrons are generated at any distance from the source as long as the energy 
of the primary cosmic rays remains in excess of the threshold for photopion 
production. This means that most gamma rays arc injected in the first few 
interaction lengths from the source, and their cascade starts at the production 
point. It is therefore clear that the result provided by the analytical calculation 
should be taken only as the estimate of an order of magnitude of the expected 
effect. Another complication consists in the lack of complete development of 
the cascade for those parts of the cascade initiated too close to the Earth. All 
these effects, included in our detailed calculations, tend to decrease the ffux of 
lower energy gamma rays compared with the rough estimate that we report 
in this section. On the other hand, while this estimate ignores the infrared 
background, some energy is transferred to the cascade due to the interaction 
with this third photon background. 

In order to have some reference numbers, we assume to have a source of UHE- 
CRs with a luminosity Lcr, in the form of cosmic rays with energy above 
some value, that we take here to be 10^^ eV. Not all of this luminosity can be 
converted into gamma rays because of the kinematic threshold in the photo- 
pion production process. If the injection spectrum is taken as E^'^'^ [4], the 
energy that can potentially be converted to gamma rays is about 20%Lcr. 

Inserting numbers in the previous expressions for the cascade spectrum, we 
obtain the following estimate of the flux in the energy region below E^^^ — 



In the TeV region, the expected fluxes are comparable with the sensitivities 
of upcoming Cerenkov gamma ray experiments such as VERITAS and HESS. 



140 TeV: 



F{E^) ^ 7.4x10 



14 




-1 
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3 The background radiation field 

In this section we describe the photon background adopted in our calculations. 
For the sake of clarity we discuss separately the cases of the radio, microwave 
and IR-Optical-UV parts of the spectrum. It is worth stressing the intrinsic 
differences between the CMBR on one side and the other backgrounds. The 
former is of cosmological origin and is pretty well known to have a Planckian 
spectrum with temperature at the present cosmic time of T = 2.728 ±0.004 K 
[21]. The radio and UV-Optical-IR backgrounds are instead due to emission 
processes in astrophysical objects and are much more poorly known obscrva- 
tionally. In particular the knowledge of the radio background in the frequency 
range below a few MHz is very poor: our Galaxy is opaque to these radio waves 
due to free-free absorption, so that we cannot measure directly the universal 
photon background at these frequencies. We are therefore forced to rely only 
upon theoretical calculations and a few observations at higher frequencies. The 
UV-Optical-IR background is also known only within a factor of a few, due 
to the difficulty in separating the extragalactic contribution from the galactic 
emission [22]. 

The UV-Optical-IR background 

In a very general way, the fluence Si, (in units of energy per unit time per unit 
surface) in photons of frequency u contributed by sources with luminosity 
function per comoving volume (j)[L,z) can be written as follows [23]: 



where fy{z, L) is the emissivity of the sources with luminosity L at red-shift 
z. 

The UV-Optical-IR background is thought to be generated by galaxies, there- 
fore we calculate here such background as a function of redshift following the 
procedure detailed in [24]. In the calculations, we adopt a cosmology with 
Ho = 71 km/s/Mpc, fi^/i^ = 0.135 with h = i/o/100 and Q = 1 (thus 
Qa = 0.73) as derived by WMAP [25]. 

The procedure described in [24] is based on building a set of synthetic galaxy 
spectral energy distributions (SEDs) as functions of the luminosity in the ISO 
band LW3 centered at 15 ^m. These spectra are built as linear combinations 
of the synthetic spectra of four galaxies, Arp 220, NGC 6090, M 82 and M 51, 
as prototypes of ultra luminous infrared, luminous infrared, starburst and 
normal galaxies respectively, as obtained from the online library of synthetic 
spectra created with the GRASIL code [26]. In this way we compute the 
emissivity fiy{z,L) at 2; = 0. The red-shifted emissivity is derived by scahng 
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the luminosity by the function g{z) of equations (12) and (14) so that: 

Mz,L) = MO,L/g{z)). (8) 



We divide the luminosity interval of the 15 /xm luminosity function of [27] in 30 
bins and for each bin we take the SED computed with the linear combination 
of the synthetic spectra. We use the best fit combination of luminosity and 
density evolution computed by [24] so that the luminosity function reads 

0(L,^)= 71(^)0(^,0) (9) 



(10) 



with 

with N = 10-°•22Mpc"^ a = 0.1, b = 0.1 [27,28]. The red-shift evolution is 
described by the quantities 

n{z)^{l + zr (11) 
giz)^il + zf^ (12) 

for z < zt and 



l + z 



0.1 



n{z)^n{zT)[^^^] (13) 



\ + z 



for z^ Zt. The parameters zt^ O-k and are as computed in [24]. 

The IR-UV spectrum has an infrared luminosity between 3.5 and 1000 /im of 
34 nWm~^sr^^. The observational bound is 50 ± 25 nWm~^sr~^ and other 
theoretical estimates lay in the range between 23 and 73 nW m~^sr~^ (see 
[29] and references therein). A direct comparison of our background with that 
obtain in [30] (dashed line) is shown in Fig. 2. 

The radio background 

As pointed out before, the universal radio background (URB) is not very well 
known mostly because it is difficult to disentangle the Galactic and extra- 
galactic components and because of the free-free absorption in the galactic 
disk and halo at low frequencies. Observations have provided us with an es- 
timate of the URB [31], while an early theoretical estimate was given in [32]. 
More recently, an attempt has been made to calculate the contribution to 
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the URB from radio-galaxies and AGNs [33] (sec Fig. 1). The issue does not 
seem to be settled, but unfortunately it seems to be confined to the field of 
theoretical calculations, at least in the frequency range below ~ 1 MHz. 




log(i/ / GHz) 

Figure 1. Figure from Ref. [33]. Contributions of normal galaxies (dotted curves), 
radio galaxies (long dashed curve), and the cosmic microwave background (short 
dashed curve) to the extragalactic radio background intensity (thick sohd curves) 
with pure luminosity evolution for all sources (upper curves), and for radio galaxies 
only (lower curves). The dotted band gives an observational estimate of the total 
extragalactic radio background intensity [31] and the dot-dash curve gives an earlier 
theoretical estimate [32]. 

In our calculation we adopt the URB as given in [33] for the case of pure 
luminosity evolution of normal galaxies and radio-galaxies. 

The local photon background obtained as a superposition of the IR-Optical- 
UV, the radio and the microwave background is plotted in fig. 2, where the 
dashed curve shows, for comparison, the result of Ref. [30] . We checked that 
our results do not change appreciably by changing the IR background, the 
reason being that most part of the cascade develops on the URB and on the 
CMB radiation. 



4 The Electromagnetic cascade: a detailed calculation 

In this section we describe our calculation of the electromagnetic cascade ini- 
tiated by a proton with ultra-high energy through photopion production off 
the photons of the CMBR. In §4.1 we describe the calculation of the rate of 
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Radiative background 




E[eV] 



Figure 2. The local photon background we used. The main peak is due to the CMB, 
the radio background is taken from [33] in the hypothesis of red-shift evolving 
sources, the IR-UV background is taken from [24]. The dashed curve shows the 
result obtained in [30]. 

injection of gamma rays and electrons in the cascade due to discrete interac- 
tions of photopion production. In §4.2 we describe in detail the propagation 
of electrons and gamma rays in the background radiation field. 



4-1 Gamma rays and electrons from photopion production 



At energies larger than the threshold for photopion production, given by 

= ' y^^) 

where e is the energy of the target photon, a proton can generate a pion in the 
final state, in the reaction of photopion production as written in Eq. 1. The 
decay of the 7r° or 7r+ generates either gamma rays {tt^ —>■ '-ff) or electrons 
(vr"'" /i^t'/i; /X — > e + z/^ + fe) . Both these secondary products can initiate a 
cascade in the universal photon background calculated in the previous section. 

The propagation of UHECRs from a source at a given distance is simulated 
using the Montecarlo approach presented in [4]. However, in the version of 
the simulation adopted here we improved the code by simulating the single 
proton-photon interactions using SOPHIA [34,35]. This event generator treats 
an interaction either via baryon resonance excitation, one particle t-channel 
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exchange (direct one-particle prodTiction), diffractive particle production and 
(non-diffractive) multiparticle production using string fragmentation. The dis- 
tribution and momenta of the final state particles are calculated from their 
branching ratios and interaction kinematics in the center of mass frame, and 
the particle energies and angles in the lab frame are calculated by Lorentz 
transformations. SOPHIA takes care also of the decay of all unstable particles 
(except for neutrons) using standard Monte Carlo methods of particle decay 
according to the available phase space. The neutron decay is implemented 
separately due to its much longer lifetime. The SOPHIA event generator has 
been shown to be in good agreement with available accelerator data. A de- 
tailed description of the code including the sampling methods, the interaction 
physics used, and the performed tests can be found in Ref. [34]. 

The gamma ray and positron yields used in our calculation are computed by 
simply accounting for the decay of neutral and charged pions. The spectra 

of gamma rays at different distances from the source, for a given distance 
between the source and the observer are plotted in Fig. 3. The spectra of elec- 
trons/positrons coming from the decay of charged pions for the same distances 
from the source are shown in Fig. 4. 

As shown in [4], the injection spectrum of UHECRs that best fits the observed 
data of AGASA and HiRes in the energy region below 10^° eV is a power law 
E~'^ with 7 = 2.6 in the case of sources that do not have luminosity evolution 
with redshift. The injection spectrum flattens to E~'^-^ for the case of sources 
with strong luminosity evolution L{z) oc (1 -|- z)^. 

In [12] a detailed analysis of the consequence of small scale anisotropics in 
the directions of arrival of UHECRs was carried out. It was found there that 
the number of doublets and triplets observed by AGASA suggests that the 
sources of UHECRs are powerful relatively rare sources, with a typical number 
density 10~® — 10~^ Mpc^'^ and a luminosity at energies above 10^^ eV of 
~ 2 X 10^2 _ 2 X 10^3 s-i_ 

The reactions of photopion production during the propagation of UHECRs 
from the sources to the Earth generate at any time new gamma rays and 
electrons, that start new cascades that will eventually reach the Earth. 

4-2 The development of the electromagnetic cascade 

As already discussed at length in section 2, the electromagnetic cascade is 
mainly driven by Inverse Compton Scattering (ICS) of high energy electrons 
on the low energy background photons and Pair Production (PP) of high 
energy photons on the low energy background photons. Higher order QED 
processes become important only for energies exceeding 10^^ eV and have a 
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18 18.5 19 19.5 20 20.5 21 21.5 

loglO(E) [eV] 

Figure 3. Photon injection spectra accumulated in bins of 10 Mpc width from the 
source (from top to bottom). 




loglO(E) LeVJ 

Figure 4. Positron injection spectra accumulated in bins of 10 Mpc width from the 
source (from top to bottom). 

significant impact only in the absence of magnetic fields and for a narrow 
peaked background spectrum [36]. In our calculations we only include ICS, 
PP and synchrotron emission in the presence of an intergalactic magnetic 
field. 

Here we use the following notation: fe{E, t) is the differential spectrum of elec- 
trons at time t during the evolution of the cascade, f^{E,t) is the differential 
spectrum of photons at the same time t. We calculate these spectra by solving 
the coupled Boltzman equations that describe the evolution of fe and due 



13 



to energy losses. 

The Boltzmann equations read [40]: 

^^^^ = -J den,{e)Rics{E,e)UE,t) + 

I demie) I dE'Pics,efe{E',t) + J denb{e) J dE'Pppf^{E' ,t) + qe{E,t) 

(16) 



dME,t) 
dt 



^-J denb{e)Rpp{E,e)f^{E,t)+ 



I deribie) J dE'PicsME', t) + qj{E, t), (17) 
where we introduced the angle averaged cross sections as 

R{E, e) = £ df } ~ f'^' <j{s{e,, E, /.)) (18) 
and the angle averaged differential cross sections as 

PiE'; E, e) = C ' fe,, da{sie,^E, „)) _ ^^^^ 
J —1 2 ci-C/ 

Here s{eb, E, fi) is the center of mass squared energy for a particle with energy 
E scattering a photon with energy at angle having cosine fi. nb(e) is the 
spectrum of background photons. qe{E, t) and q^{E, t) are the source functions 
for electrons and photons respectively, calculated as described in §4.1. At any 
time new electrons and gamma rays are generated and start new cascades, 
due to the reactions of photopion production of protons. 

4-2.1 Pair production 

The pair production cross section is given by [37]: 



<7pp = (l - P') [(3 - P') In 1±| - 2/? (2 - f3) 



(20) 



where P is the velocity of the outgoing electron in the Center of Mass (CM) 
frame. 

The threshold for pair production is: 

"^1; o c mil ^ ^ ^ 



Etn^^^ 2.6 X 10^^ ( — ) eV. (21) 
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The spectrum of electrons generated by pair production in the interaction of 
two photons with energy and with E^ ^ is [38]: 



dN^.^ 3 1 



64e2£;3 [{A 



A A ( E'\ 



2{2eA-l) 



(A-E',)E', 



+ 



1 7 



A^ 



eAj (A - E'.YEf 



(22) 



where all energies are in units of the electron mass and we put A = e + E^. 
The range of electron energies in the final state is defined by: 



A 
2 



1- 1 1- 



1 



A 



1 + Wi- 



se 



(23) 



^.2.2 Inverse Compton Scattering 

The cross section for ICS {cfb — > 67) is given by the well-known Klein-Nishina 
expression (e.g. [39,40] ) : 



3 ml 

cries = ^T--—s 
8 sp 



r^(2 + 2/3 - - 2/3^) - j^(2 - 3/3^ - /3^) In 



/3(l + /3)^ ' " " " ^ /32^ " " ^ 1-/3 

(24) 

where /3 = (s — m^) /(s -|- mg) is the velocity of the outgoing electron in the 
center of mass frame. 

As in the case of pair production, we use an approximate formula for the 
spectrum of the produced photons/electrons. The approximation holds as long 
as the incoming electron is relativistic. The spectrum of electrons resulting 
from the scattering is [41,42]: 



dA^^ / ^ ^, X 3 



2q + logg + (1 + 2q){l -q) + \ {Teqf 



(25) 

where 7 = ^1 + El, Ei = E'^/j, Tg = 4e7, q = p^^^^^^ and the range is 

restricted to J- < £^1 < Y+p-- All energies are again in units of the electron 
mass. 

The differential spectrum of the outgoing photons is obtained by substituting 
E^-E'^ to ineq. (25). 

By averaging the spectra of equations (22) and (25) over the background 
spectrum we can obtain directly the integrated quantities of Eq. (17): 

/r (IN 
denie)PppiE'^; E„ e) ^ j den{e)^ {E„ e, E'J (26) 



15 



and 

den{e)PicsAE',; Ee, e) ^ J den(e)^ (e^, e, E'^) (27) 
and similarly for the outgoing electrons. 



4 -2. 3 The magnetic field 

Synchrotron losses of high energy electrons are introduced as continuous en- 
ergy losses. This works perfectly for the range of energies that we are interested 
in. 

The rate of energy losses can be written as: 

dE^ 4 f E \\ . ^ 17/ . , f f E \^ GeV 

— ^ = — car — ^ (sin a) ~ 7.9x 10"^^ (sin a — — — 

dt 3 Stt Vmec2; ^ ' ^ ' VnGy V GeV J year 

(28) 

where a is a pitch angle. 

The spectrum of photons radiated by synchrotron emission by an electron 
with energy E, assuming isotropy of the electron distribution, is [39] 

dP 2y/3 ^ ( E'\ 

where — 9.27 x 10~^^ erg/G is the Bohr magneton, a is the fine structure 
constant, 

= 3/iB572 ~ 6.6 X 10-20-^ GeV (30) 

is the critical energy and 

POO 

F{x)=x K,/,{i)di (31) 

Jx 

with K^i^{^) the modified Bessel function of fractional order. The power spec- 
trum peaks at ~ 0.29£'c. 



5 Results 



In this section we discuss the results of our numerical calculations of the devel- 
opment of the electromagnetic cascade from nearby sources of UHECRs. We 
first discuss the case in which the cascade develops in the absence of magnetic 
field in the intergalactic medium (§5.1). In §5.2 we apply our calculations to 
a magnetized universe. 
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Figure 5. Flux of gamma ray at the Earth as resulting from the injection of cosmic 
rays in continuous point sources located at the distances indicated in the graph. 
The thick hues are the sensitivities of MAGIC, HESS, VERITAS and GLAST. The 
left (right) panel is for a source with luminosity 2 x 10^'^ ergs~^ (2 x 10*^^ ergs^^) 
at energies larger than 10^^ eV. 



5.1 The unmagnetized case 



Numerical calculations of the propagation of UHECRs from point sources ap- 
pear to give small scale anisotropies compatible with the observations carried 
out with the AGASA experiment for a source density of 10~^ — 10~^ Mpc~^, 
corresponding to a source luminosity above 10^^ eV of 2 x 10^^ — 2 x lO^^erg s~^. 
The best fit injection spectrum is E~'^'^ for sources with no luminosity evolu- 
tion and for sources with strong luminosity evolution, L{z) oc (1 + z)'^. 
We followed the propagation of UHECRs with the simulation used in [4,12] in 
which the interactions are implemented with the SOPHIA package. At each 
step of the simulation, the spectra of gamma rays (from the decay of neu- 
tral pions) and of electrons (from the decay of charged pions) are determined 
and used as source functions for the gamma rays and electrons respectively in 
the coupled Boltzmann equations for the development of the electromagnetic 
cascade. 

The results of our calculations for a proton injection spectrum are shown 
in Fig. 5 for a source located at 6 Mpc (short dashed line), 21 Mpc (solid line), 
48 Mpc (long dashed line) and 71 Mpc (dash-dotted line). In the same plot we 
show the sensitivities of MAGIC, HESS, VERITAS and GLAST (thick lines). 

For the case with higher luminosity (left panel) the predicted fluxes should 
be detectable by the experiments HESS and VERITAS for all the considered 
distances. MAGIC should detect a signal only if the source is at about 20 
Mpc distance. The corresponding fluxes for the case with source luminosity 
2 X 10^^ ergs~^ are slightly below the sensitivities of these telescopes. None 
of the predicted fluxes appears to be detectable in the lower energy range. 
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Figure 6. Flux of gamma ray at the Earth as resulting from the injection of cosmic 
rays in continuous point sources located at the distances indicated in the graph. 
The thick hues are the sensitivities of MAGIC, HESS, VERITAS and GLAST. The 
left (right) panel is for a source with luminosity 2 x 10^^ ergs~^ (2 x 10*^^ ergs^^) 
at energies larger than 10^^ eV. The luminosity evolution of the sources is taken as 
L{z) oc (1 + z)^ and the injection spectrum of UHECRs as E''^-^. 



accessible to GLAST. It is worth to remind that the largest distance considered 
in the calculations is comparable with the loss length of UHECRs with energy 
in excess of 10^° eV, for which a counterpart is still lacking. 

The flux of gamma rays does not necessarily decrease for larger distances to the 
source. This is due to the fact that for too small distances, the electromagnetic 
cascade has not enough time to develop and the lower energy part of the 
spectrum does not reach the equilibrium shape. This is the reason why in Fig. 
5 the spectrum of gamma rays from a source at distance 6 Mpc (short dashed 
line) for a wide range of energies is below that corresponding to a source at 
larger distances. 

Another point to keep in mind is that the calculation in [12] were carried out 
for a class of sources with identical luminosity. It would probably be more 
realistic to imagine that the sources of UHECRs have some kind of luminosity 
function, and that the value used here is some sort of average luminosity. In this 
case it is quite plausible that sources with cosmic ray luminosity larger than 
2 X 10^^ erg s~^ exist in nature. Since the flux of gamma rays in the cascade 
scales linearly with the cosmic ray luminosity, the cascade signal increases 
correspondingly. 

For sources with a strong luminosity evolution, the flatter required injection 
spectrum in the form of UHECRs slightly increases the predicted cascade 
fluxes, as shown in Fig. 6. 
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5.2 The case of a magnetized universe 

The development of electromagnetic cascades in the universe is very sensitive 
to the presence of magnetic fields, mainly in two ways: 1) the presence of a 
magnetic field modifies the duration of the cascading processes as observed 
at the Earth, due to the deflection of the charged components in the cascade; 
2) electrons with very high energy can lose energy by synchrotron emission, 
therefore subtracting energy from the cascade. 

The first effect is important only for the case of bursting sources. In that 
case the duration of the burst in the form of the electromagnetic cascade may 
depend on energy, namely the duration of the burst appears to be different at 
different energies. Very tight limits on the magnetic field in the IGM can be 
imposed in principle from the detection or non detection of the lowest energy 
part of the cascade [43]. Since we consider here continuous sources, we will 
not investigate any further this first effect of the magnetic field, and we will 
concentrate our attention on the role of synchrotron losses of the electron 
component. 

At very high energies, the cross section for ICS is in the Klein-Nishina regime. 
This is the energy region where it is easier for the synchrotron losses to dom- 
inate over ICS. Therefore we expect that the role of the magnetic field is 
particularly relevant when the field is concentrated around the source, where 
the electrons in the cascade have the highest energies. Moreover, as discussed 
in §4.2.3, synchrotron losses may become dominant on ICS in the energy region 
around 10^° eV if the field is larger than ~ 10~^° G. 

Our results for the case of a source with luminosity 2 x 10^^ erg s^^ are 
summarized in Fig. 7 for four values of the magnetic field, namely 10~^^ G, 
10~^° G, 10~^ G and 10~*^ G. The distances to the source are taken as in Fig. 
5, while the source luminosity in the form of cosmic rays with energy above 
10^^ eV has been taken to be 2 x 10"^^ erg s~^. 

The effect of the magnetic field on the development of the electromagnetic 
cascade is clear: for magnetic fields smaller than 10~^ G the predicted fiuxes 
start to be at the level of sensitivity of VERITAS and HESS, and are unde- 
tectable for larger magnetic field strengths. On the other hand, one has to 
keep in mind that a magnetic field of 10~^ G widespread in the whole IGM 
is already the highest possible value compatible with Faraday rotation mea- 
surements [44]. Larger fields are possible only on more compact regions: for 
this reason we also explore the situation in which a magnetic field is present 
only in a neighborhood of the source or near the Earth. A recent investigation 
of the effects of this magnetic field on the propagation of UHECRs has been 
presented in [45,46]. In particular in Ref. [45] it was shown that the average 
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L = 2xl0" erg/s, B = 10"' ' G L = 2x10''"' erg/s, G = 10"'° G 




Energy [GeV] Energy [GeV] 

Figure 7. Flux of gamma rays at the Earth for a magnetic field in the inter- 
galactic medium B = W~^^, 10-^°, IQ-^, 10"'^ Gauss and a source with luminosity 
2 X 10^'^ ergs~^ at energies larger than 10^^ eV. 

magnetic fields are quite smaller that 10^^ G, with the exception of the fila- 
mentary regions where they are of the order of 10^^ G. Clearly much higher 
fields are achieved in Mpc scale virialized regions such as clusters of galaxies. 
This conclusion does not appear to be confirmed in [46]. 

In the left panel of Fig. 8 we plot the fluxes of gamma rays calculated if the 
source is located in a region of 10 Mpc where the magnetic field is 10^^ G. 

Similarly, in the right panel of Fig. 8 we plot the fluxes of gamma rays calcu- 
lated in the case that the Earth is located in a region of 10 Mpc in which the 
field is 10^^ G. This might be the case for the Local Supercluster. 

It is easy to see that if the magnetic field is concentrated close to the detector 
and not close to the source, its effect on the development of the cascade is 
negligible, as can be inferred by comparing the right panel of Fig. 8 and the 
left panel of Fig. 5. 

In the left (right) panel of Fig. 9 we plot the expected fluxes in the case that the 
magnetic field is 10~^ G in a region of size 1 Mpc around the source (Earth). 



20 




Figure 8. Flux of gamma rays at the Earth for a magnetic field in the intergalactic 
medium B = 10^^ Gauss located in a region of 10 Mpc around the source and the 
Earth respectively. The source luminosity is 2 x 10'*'^ ergs~^ at energies larger than 
IQi^ eV. 



L = 2x10 erg/s L = 2x10 erg/s 




Figure 9. Flux of gamma rays at the Earth for a magnetic field in the intergalactic 
medium B = 10^^ Gauss located in a region of 1 Mpc around the source and around 
the Earth respectively. The source luminosity is 2 x 10^^ ergs~^ at energies larger 
than 10^9 eV. 

The expected fluxes in these cases are not much affected by the presence of 
the magnetic field: in the near source region of 1 Mpc almost no reaction 
of photopion production occurred, therefore no cascade was developed yet. 
In the near Earth region the cascade has already reached in most cases the 
final stages, where the effects of synchrotron losses of electrons are irrelevant 
compared with ICS. 



6 Conclusions 



The sources of UHECRs are still unknown mainly because we do not know 
yet whether they are generated in a few very powerful sources or rather in 
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numerous low power sources. This degeneracy has probably been broken after 
the discovery of small scale anisotropics in the arrival directions of UHECRs in 
the AGASA experiment, although the statistical significance of this discovery 
appears to be still controversial [10]. In [12], the AGASA data were used to 
infer that the most likely number density of sources of UHECRs is around 
10~^ — 10~^ Mpc~^, which corresponds to a source luminosity of 2 x 10^^ — 
2 X IC^^ erg s^^ at energies above 10"'^'^ cV. These estimates, obtained in the 
absence of intergalactic magnetic fields, remain valid as long as the magnetic 
field does not exceed ~ 10~^° G over cosmological distances. 

During their propagation from the sources to the Earth UHECRs suffer pho- 
topion production if their energy is above threshold for this process, therefore 
generating charged and neutral pions. While the former decay mainly into elec- 
trons and neutrinos, the latter generate gamma rays. Both gamma rays and 
electrons initiate an electromagnetic cascade that transfers energy from high 
energies to lower energies, more easily accessible to gamma ray telescopes. We 
calculated the development of this electromagnetic cascade with and without 
the presence of an intergalactic magnetic field, and investigated the detectabil- 
ity of the gamma ray signal with ground based telescopes such as MAGIC, 
VERITAS and HESS, as well as with space-borne telescopes such as GLAST 

4 

The cascade radiation appears to be detectable by VERITAS and HESS for a 
source luminosity around 2 x 10^^ erg in the form of cosmic rays with energy 
larger than 10^^ eV, if the magnetic field in the IGM is constant over all the 
propagation volume and smaller than 10~^ G (This luminosity corresponds to 
the lower source density compatible with the small scale anisotropics observed 
by AGASA). If the magnetic field is patched, stronger fields are allowed close 
to the Earth without appreciably changing the previous conclusions. If the 
region with stronger fields is close to the source location, the initial stages 
of the cascade may be affected if the magnetized region is large enough, so 
that the energy content of the cascade gets suppressed. A large field close 
to the Earth does not affect the conclusions in any appreciable way, since 
the cascade is already fully developed at that point. A magnetic field of InG 
near the source implies slightly lower gamma ray fluxes if it is extended in 
a region of about 10 Mpc. A stronger fleld, of say 10 nC, in a 1 Mpc region 
either around the source or the observer does not change the results in any 
appreciable way. 

The main uncertainty involved in our calculations is the luminosity of the 
sources in the form of UHECRs. The numbers we adopted are derived in [12], 
on the basis of the appearance of small scale anisotropics in the AGASA data 



^ The fluxes derived in the present paper could be higher by up to a factor of ~ 2 
due to the contribution from proton pair production, that was neglected here. 
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and in the assumption of sources that accelerate UHECRs continuously. For 
bursting sources, the luminosity per source may be many orders of magnitude 
larger than that adopted here, but it would be concentrated within a short time 
interval. The diffuse and time-spread appearance of the observed UHECRs is 
achieved, within this class of models, by assuming the presence of a small 
magnetic field in the intergalactic medium, responsible for sizeable time delay 
that make the burst look like a continuous signal when observed through 
charged particles. This is the case for the gamma ray burst models [47,48,49] 
of UHECRs. 

A similar time delay is introduced in the development of the electromagnetic 
cascade. Moreover, an additional time delay is introduced due to the presence 
of electrons and positrons in the cascade itself. These effects cannot be ac- 
counted for within the framework proposed here, and will be considered in a 
forthcoming publication. 
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